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We have studied the ordering of a two dimensional charge 
stabilised colloidal system in the presence of a stationary one 
dimensionally modulated laser field formed by the superposi- 
tion of two modulations with wavevectors q^T and qo/r, where 
go is the wavevector corresponding to the first peak of the di- 
rect correlation function of the unperturbed liquid, and r is 
the golden mean. In the framework of the Landau- Alexander- 
McTague theory we find that a decagonal quasicrystalline 
phase is stabler than the liquid or the triangular lattice in 
certain regions of the phase diagram. Our study also shows a 
reentrant melting phenomena for larger laser field strengths. 
We find that the transition from the modulated liquid to the 
quasicrystalline phase is continuous in contrast to a first or- 
der transition from the modulated liquid to the triangular 
crystalline phase. 

PACS numbers: 82.70.Dd, 64.70.Dv 

In recent years, there has been considerable interest in 
the ordering of charge stabilised colloidal particles in the 
presence of stationary laser modulations Due to 

their large diameter and charge (^lOOOe for a particle 
of diameter 1000 A), charge stabilised colloidal particles 
have large polarisabilities. The electric field of the laser 
beam induces dipole moments on the colloidal particles 
and these dipole moments in turn interact with the laser 
electric field. The resulting interaction energy of a col- 
loidal particle at position r is equivalent to an external 
potential Ve{r) = —^x{E{^))^j where % is the dielectric 
susceptibility of the colloidal particles and E{r) is the 
electric field at r. This interaction causes the particles 
to preferentially sit at the maxima and the minima of 
the electric field modulations and thus promotes density 
modulations at twice the wave-vector of the modulating 
electric field. By manipulating the field pattern one can 
get complex structures, referred to as Optical matter [0, 
for moderate field strengths which are easily attainable 
in the laboratory (for recent reviews on colloids and other 
ref.s see When the interaction energy with the exter- 
nal field is of the same order as the thermal energy of the 
particles, ordering can take place with wavevectors other 
than the wavevectors of the periodic stationary external 
potential Ve{r). This is due to the nonlinear coupling 
of the different order parameter modes. In particular, it 
has been observed that when the wavevector of an one 
dimensional laser modulation is tuned to the wavevector 
go corresponding to the first peak of the direct correlation 
function of the unperturbed liquid, the system freezes to 
a triangular lattice and simulation studies have 
also reproduced the experimental findings. Several is- 



sues connected with this phenomenon of " Laser Induced 
Freezing (LIF)" have been explored in a recent Density 
functional theory Jp|. 

Burns et al. [|l| have demonstrated that a quasi- 
crystalline arrangement of colloidal particles in two di- 
mensionally confined geometry is obtainable by subject- 
ing them to a superposition of five equiangular coherent 
laser beams. In their experiment the external field had 
the same symmetry of the final structure. A two dimen- 
sional quasicrystal is characterised by density modula- 
tions in four independent wavevectors It is interest- 
ing to ask whether one can generate LIF into a quasicrys- 
talline structure: i.e., have the external modulation cou- 
ple directly only to a few of the order parameter modes, 
and have the system generate density modulations with 
the other wavevectors needed to make up the quasicrys- 
talline structure via the nonlinear coupling of the order 
parameter modes. 

In this work, we show using the Landau-Alexander- 
McTague theory that this is indeed possible, using 
a superposition of two 1-d laser modulations with their 
wavevectors tuned to be qqt and qo/r, where r(= '^2^ ) 
is the golden mean. We show that in the presence of 
such an external modulations the system undergoes a 
transition to a quasicrystalline structure having decago- 
nal symmetry for a certain range of parameters. We also 
observe a reentrant liquid phase (see figure ||). 

The Landau- Alexander-McTague theory is a mean 
field theory which expresses the free energy(jr) as a poly- 
nomial expansion in powers of the order parameters(pi) 
of the system; where pi are the fourier components of the 
density, (or better, of the molecular field) : 

+D[Y^ pjf + E[Y^ pi] - ye{G,)pr-^ (1) 

i i i 

Here all the parameters Bi,C,D and E are assumed 
to be positive, and Ve{Gi) are fourier coefficients of the 
external modulation. To find the most stable configura- 
tion of the system for a given set of parameter values, 
one minimises this free energy for different choices of the 
order parameter sets corresponding to different lattice 
arrangements. The quadratic term in the free energy 
favours the liquid phase and the quartic terms ensure 
global stability. When Bi are sufficiently small, the cu- 
bic term clearly drives the system towards structures in 
which several sets of three wavevectors obey a triangular 
relationship among themselves. In the absence of exter- 
nal modulations, i.e., of Ve{Gi), in 2-d this would favour 
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a triangular lattice. A decagonal quasicrystal would not 
be a stable phase, because no triangular relationships ex- 
ist among the wavevectors corresponding to a decagonal 
symmetry. 
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G„ = tG, 
Of, = (-l/T)Gi 
-c = (5" + l)/2 



FIG. 1. Density wave- vectors corresponding to (a) modu- 
lated quasicrystal and (b) modulated triangular lattice. Note 
that for each density fluctuation at wavevector G, density 
fluctuation at - G also is present, though not shown in the 
figure. 

Now consider what happens when an 1-d quasiperi- 
odic external laser modulation potential of strength Ve 
at wave vectors ±Gq, and ±G/3 , where Ga — tGi and 
G/3 = — — , and Gi is any wavevector with magnitude goi 
is applied. Our choice of wavevectors is motivated from 
the fact that these wavevectors have triangular relation- 
ships with several wavevectors Gi of magnitude go and 
corresponding to a decagonal symmetry [Figure |l|.(a)]. 
Noting that r — = 1 , we get the following three triangu- 
lar relations among the wave vectors defining a decagonal 
quasicrystal. 



G^^ — Gq^ — G — 
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By contrast, for a triangular lattice one finds two such 
sets of wave- vectors which add to zero [figure 0.(b)] in 
the presence of the same external field: 



G^ ^ G(x — 0? ^1 
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g(*)-g^*) = o 
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Note that all of the triangular relations for the decago- 
nal quasicrystal necessarily involve the wavevectors of the 
external modulation potential. In contrast, only one of 
the two relations for the triangular lattice involves the 
modulation wavevectors. 

Using the above relations, we can write down the free 
energy for a decagonal quasicrystal in the presence of our 
modulating potential as: 
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Here, using symmetry 
pf=pf ^^^pf^pf 



considerations we have set 
Thus we now have five inde- 
pendent order parameters for the modulated decagonal 
structure. Accordingly we get five coupled polynomial 
equations from minimising the free energy with respect 
to these order parameters. These we have solved numeri- 
cally using Newton's method of finding roots, for varying 
values of the parameters i?o and and fixed values, 
B\ = 0.15, G = 1.0, D = 0.125 and E = 0.75, for the re- 
maining parameters. The fixed parameter values chosen 
here are the same as in ||^ ; except for Bi , which we have 
chosen for convenience to be the same at Ga and G/3, and 
to be much larger than the maximum value considered 
for Bq. This is motivated from the fact that the liquid 
structure factor peaks at go and having a density modu- 
lation at some other wave-vector will cost more energy. 
The qualitative features of our results are not sensitive 
to these specific values. 

A similar, but separate, calculation is done for the tri- 
angular lattice with two independent order parameters 
corresponding to triangular symmetry, (since, by sym- 
metry considerations, ^2*"* = Ps*"* ) and the remaining two 
order parameters corresponding to the external modulat- 
ing field. The corresponding free energy is given by: 

T^^^^-2VM^+p^^]+2Bo[pf+2pf] + 

2Bi[pr + pf ] - 2G[prVr^V + 
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FIG. 2. Free energy difference for Bo = 0.02. The free 
energy which is lower than others corresponds to the stable 
phase at a particular external field. 

In Newton's method for root finding, when multiple 
solutions, corresponding to multiple minima (more gen- 
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erally extrema) of the free energy, are present, the so- 
lution to which the result converges depends on the ini- 
tial inputs for the order parameters. In our calculations 
each of the two free energy functions corresponding to 
the triangular lattice and the quasicrystalline structure 
was effectively minimised with the initial guess value for 
the induced order parameters being 0.0 or 0.5. We found 
that the results converged to one of the two solutions: one 
corresponding to the modulated liquid phase, where only 
Pa{= Pp) and pi are nonzero, and the other to the modu- 
lated crystalline or the quasicrystalline phase, where the 
other order parameters are also nonzero. Typically, a 
maximum 20000 iterations were carried out in Newton's 
method for finding roots. It was checked that the results 
were not dependent on the initial guess values except for 
the flow to one of the two solutions alluded to above. For 
each parameter set the phase corresponding to the lowest 
free energy amongst the different solutions was chosen as 
the stable phase. 
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FIG. 4. Order parameters corresponding to triangular or- 
der (p'*' = pj*' = Pa*') a^rid those corresponding to quasicrys- 
talline order (p(«' = p^"^ = p^"' = pf'> = p^"') for Bo = 0.02. 
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FIG. 3. Results for the order parameters (p^ = p^), di- 
rectly coupling to the external field and for the order param- 
eter pi with wavevector Gi parallel to the wavevector of the 
external modulation for Bo = 0.02. 

In Figure ^ we present results for the free energy differ- 
ences for the three different kinds of order as a function of 
Ve for Bo = 0.02. From these the phase with the lowest 
free energy is easily chosen, and the corresponding order 
parameters are shown in figures ^ and ^ as functions of 

The full phase diagram is shown in figure |^ . We term 
the phases as modulated, since the order parameters with 
wavevectors along the wavevectors of the external mod- 
ulation are typically higher than those with wavevectors 
in the other directions. 

Depending upon the value of , as the external field 
is increased from zero, several interesting transitions are 
discernible from the phase diagram. For low values of Bq 
and of the modulation potential, the system is a triangu- 
lar lattice. As the field strength is increased, the system 
undergoes a first order phase transition to the modulated 



decagonal quasicrystal. At still higher field values it goes 
to the modulated liquid phase via a continuous transition 




FIG. 5. Phase diagram from the Landau - Alexander - Mc- 
Tague theory 

For i?o between 0.02 and 0.025, with increasing field 
strength one encounters a first order LIF transition from 
the modulated liquid to the triangular crystal; a first or- 
der transition from the triangular crystal to the decago- 
nal quasicrystal and finally a continuous transition to 
the modulated liquid reentrant phase from the quasicrys- 
talline phase (fig. 

It is worth pointing out that near — 0.025 one can 
go from the modulated liquid directly to the quasicrys- 
talline phase via a continuous LIF transition and thence 
to the triangular phase via a first order transition. With 
increasing field the system goes once again to the qua- 
sicrystalline phase through a reentrant first order transi- 
tion. At still higher fields the transition to a reentrant 
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liquid phase occures via a second order phase boundary. 

For Bq > 0.0265 the system goes into the quasicrys- 
talUne order from the hquid phase via a continuous tran- 
sition as the field strength is increased and then again 
melts to the modulated liquid phase at still higher field. 
Above Bq — 0.05 the system remains in the modulated 
liquid phase for all field strengths. 

The fact that at low external fields, the system either 
freezes to a triangular lattice or remains a liquid depend- 
ing upon the value of Bq and that the transition is first 
order is consistent with our knowledge about colloidal 
systems in the absence of any external field modulation. 

At moderate field strengths, the system gains energet- 
ically by having density modulations corresponding to 
the wave- vectors of the modulating laser field. At higher 
field strengths, even though the external field modulation 
is 1-dimensional, density modulations corresponding to a 
2-dimensional triangular or decagonal symmetry develop 
because of the nonlinear coupling among the order pa- 
rameter modes. 

At still higher fields the continuous melting is con- 
sistent with previous simulations Q and the Landau- 
Alexander-McTague meanfield analysis of the laser in- 
duced freezing P,p^. 

The boundary between the triangular lattice and the 
quasicrystalline phase is first order because the two struc- 
tures are of completely different symmetry and one can 
not deform a triangular lattice continuously to get a 
decagonal symmetry. 

In contrast, the phase boundary between the mod- 
ulated liquid and the crystalline and quasicrystalline 
phases can be first order or continuous. The mechanism 
that determines which has been discussed in detail in ref 
1^ . Basically, in the modulated liquid phase pa — pp and 
pi are nonzero. Now consider setting up a Landau expan- 
sion for the free energy in powers only of the additional 
order parameters that characterise the crystalline or qua- 
sicrystalline phases ( i.e., with respect to the modulated 
liquid phase) . Such an expansion has only even order in- 
variants in the additional order parameters for the cases 
we are discussing. The order of the transition depends 
on the signs of T2 and T4, the second order and fourth or- 
der coefficients respectively, in such an expansion. A first 
order transition ensues when T4 < 0; whereas a continu- 
ous transition results when T4 > 0, the phase boundary 
being determined by the condition T2 = 0. 

When the field strength is very large, the linear term 
coupling the laser modulation potential with the order 
parameters at the modulation wavevectors seems to be 
the most dominant term for lowering the free energy.And 
we find that having density modulations along other 
wavevectors no longer lowers the free energy. So the 
decagonal phase melts to give the reentrant modulated 
liquid [|lO|. 

In conclusion we have shown that ,by subjecting a 2- 
dimensionally confined charge-stabilised colloidal liquid 
to a superposition of two 1-d laser modulations with their 
wavevectors tuned to be go''' and q^/T (where go is the 



wavevector of first peak of the liquid structure factor), 
one can generate laser-induced-freezing into a decagonal 
quasicrystalline order. We have also shown that for larger 
laser field strengths, within mean field theory this transi- 
tion is continuous and shows a reentrant melting back to 
the modulated liquid phase. It would be of interest if the 
experiments of ref. |3| can be extended to explore these 
transitions and the resulting quasicrystalline phase. 
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